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Abstract.  The  rate  of  reaction  of  ground  state 
Op  with  has  been  measured  in  the  energy  range, 
0-05-11  eV.  The  measurements  show  that  the 
r^e  coefficient  does  not  decrease  drastically 
with  increasing  energy,  contrary  to  previous  work. 
The  rate  constant  measured  here  is  then  applied  to 
a  simple  model,  which  includes  reactions  of 
with  water  dimers,  to  calculate  water  concentra¬ 
tion  in  the  vicinity  of  the  Space  Shuttle.  The 
results  suggest  that  the  pressure  of  water  vapor 
may  be  lower  than  previously  thoiwht^with 


for  the  formation  of  water  dimer  by  evaporation 
would  lead  to  better  agreement  between  the  data 
derived  from  the  neutral  mass  spectrometer 
measurements  2Lnd  those  derived  from  ion  mass 
spectra. 

Experimental  and  Results 


upper  limit  of  1  x  lol  molecule 


Introduction 


The  flight  of  the  space  shuttle  and  the  'i 

availability  of  composition  data  [Hunten  and  ' 

Calo,  1 985 :Wulf  and  von  Zahn,  1986;  Grebowsky  et 
al . ,  1987a:  Grebowsky  et  al . ,  1987b]  has  renewed 
interest  in  the  rates  of  ion-neutral  reactions  at 
high  kinetic  energies.  The  high  orbital  velocity 
of  the  space  shuttle  (7.7±  0.2  km  s"')  leads  to 
reactions  between  the  ambient  atmospheric  ions 
(principally  0*)  and  the  local  gaseous  constit¬ 
uents  of  the  space  shuttle  at  collision  energies 
of  <  4.5  eV  for  0*  (depending  on  the  neutral 
target  and  on  the  orientation  of  the  space  shuttle 
with  respect  to  the  velocity  vector). 

Understanding  the  origin  of  large  amounts  of 
H,0*  ions  in  the  bay  of  the  space  shuttle  [Naroisi 
et  al .  ,  1983]  has  not  been  straight-forward:  /T"' 

relating  the  Intensity  of  HjO'  to  the  local  •  /  , 
pressure  of  H  O  is  not  easy,  although  an  estiraatel  | 
of  [HjO]  =  0,6-2. 4  X  10’  cm"’  has  been  made  I  \ 

[Caledonia  et  al .  ,  1987].  The  complicating  V 

circumstances  are:  (1)  the  uncertainty  in  the 
interaction  energy  for  the  reaction  of  0*  with 
H  O;  and  (2)  whether  the  measured  0*  and  HjO* 
abundances  can  be  derived  from  the  mass  spectro¬ 
meter  data  since  the  velocities  (and  hence 
collection  efficiency  of  the  mass  spectrometer)  of 
0*  and  HjO*  are  likely  to  be  quite  different 

[Hunton  and  Calo,  1985,  Grebowsky  et  al . ,  1987b]. 

A  previous  measurement  of  the  cross  section  for 
the  reaction 

0*  HjO  HjO*  0  (1) 

as  a  function  of  collision  enerw  in  a  beam-target 
experiment  [Murad  and  Lai,  1986]  lead  to  a  ques¬ 
tion  as  to  the  origin  of  H.O* • in  view  of  the  large 
amounts  of  accompanying  H^O* .  By  comparison  with 
[HjO]  derived  from  ion  abundances,  results  from 
the  neutral  mass  spectrometer  [Wulf  and  von  Zahn, 
1986]  suggest  that  [H^O]  -2x10®  cm"’.  This 
work  was  undertaken  in  order  to  remeasure  the 
cross  section  for  reaction  (l)  as  a  function  of 
energy  under  conditions  where  the  state  of  0*  is 
much  better  controlled  and  the  collision  energy  is 
better  determined.  It  was  hoped  that  the  new 
measurement  and  the  introduction  of  a  mechanism 
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The  results  reported  here  were  obtained  in  a 
r  triple  cell  Ion  Cyclotron  Resonance  mass  spectro- 

lUhL^with  li^^^eter,  whose  operation  has  been  described 

^7 A  '  .jA-eviously  [Heninger  et  al.,1986:  Mauclaire  et 
al.,  1987].  In  the  present  experiment  the  three 
aTL,. cascaded,  differentially-pumped,  ICR  cells  are 

)  used  sequentially  for  the  following  purpose;  l.ion 
production:  2.0*  selection;  and  3.0*  acceleration, 
\  reaction  and  detection. 

Production  of  ground  state  0*  was  made  by  30 
eV  electron  impact  on  COj,  known  to  give  0*  (*S) 
with  negligible  kinetic  energy  [Murad,  1973]. 

.  These  primary  ions  were  selected  by  removing  all 

'  unwanted  ions  using  selective  ejections  at  cyclo¬ 

tron  frequencies  in  cells  H  1  and  H  2.  The  0*  ions 
were  drifted  and  then  stored  in  cell  H  3  for  a 
fixed  50  ms  duration.  At  the  beginning  of  this 
storage  time,  the  0*  ions  were  given  the  desired 
kinetic  energy  by  excitation  of  their  cyclotron 
motion  by  a  resonant  radio-frequency  field  applied 
-  perpendicularly  to  the  magnetic  field.  This  energy 

isi  given  by  KE  *  q®E^t^/8m,  where  q  and  m  are 

respectively  the  charge  and  mass  of  the  ion,  while 
E  and  t  the  amplitude  of  the  rf  field  and  its 
1  \duration,  respectively,  t  was  kept  fixed  at  1 
]  1  jroillisecond,  and  the  amplitude  of  the  rf  field  was 
I  \  /adjusted  to  obtain  the  desired  kinetic  energies. 

The  end  of  the  accelerating  pulse  defines  the 
origin  of  the  reaction  time,  which,  using  an  H  O 
pressure  of  1x10*^  torr,  was  varied  from  4  to  30 
ms.  and  was  followed  by  detection  of  the  ions. 
Because  of  the  fast  secondary  reactions  of  H,0* 
qj.  with  HjO,  for  each  0*  kinetic  energy  reported,  the 
reaction  rates  were  obtained  from  the  decay  of  0* 
as  a  function  of  time. 


Although  the  width  of  the  energy  distribution 
is  rather  difficult  to  ascertain,  the  validity  of 
this  experimental  procedure  has  been  tested  using 
a  reaction  whose  branching  ratio  is  strongly 
dependent  on  kinetic  energy; 


N*  0, 


0,*  N 
NO*  0 
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Our  measured  k^^/k^^^  varied  with  kinetic  energy  in 
a  manner  which  agreed  with  the  flow-drift  tube 
determination  [Howorka  et  al . ,  1980]. 

The  pressure  of  H.O  was  measured  with  an 
ionization  gauge  calibrated  by  the  determination 
of  the  rate  of  reaction  of  Kr*  with  H,0  at  thermal 
energy  [Marx  et  al.,  1983].  The  validity  of  this 
procedure  is  supported  by  the  agreement  between 
our  thermal  rate  constant  and  previous  deter¬ 
minations  . 

Figure  1  shows  the  results  obtained  here  for 
the  cross  section  of  reaction  (1)  as  a  function  of 
energy  in  the  range  »  0.035-11.5  eV.  This 
cross  section  curve  is  different  from  that  _ 
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Figure  1-Rate  constant  for  reaction  (1)  as  a 
function  of  kinetic  energy  of  0* . 


reported  earlier  [Murad  and  Lai.  1985]:  At  the 
lowest  energy  (E^»0.2  eV)  reported  by  Murad  and 
Lai.  the  ICR  cross  sections  are  nearly  a  factor  of 
3  higher  than  the  beam  experiment,  while  at  E  • 
10  eV.  the  ICR  data  is  about  a  factor  of  2  hl^er. 
The  reason  for  this  discrepancy  is  not  clear  at 
this  time,  although  it  is  possible  that  both  the 
energy  calibration  (measuring  the  point  of  maximum 
slope  in  the  retarding  potential  curve  for  the 
primary  ions)  and  the  cross  section  calibration 
(which  was  done  by  comparing  the  measured  cross 
section  with  that  of  N^D'+D)  may  be 

incorrect.  In  any  case,  we  believe  the  present 
results  are  more  reliable,  especially  since  now 
the  rate  constant  does  not  vary  much  in  going  from 
thermal  energies  to  11  eV,  as  given  in  Table  1. 


TABLE  1 


Energy.  Ej,  (eV) 

0.0327 

2.5 

230 

0.23 

2.9 

125 

0.46 

2.8 

85.5 

0.58 

2.8 

77.7 

0.92 

2.6 

57.6 

1  .46 

2.5 

42.8 

1  .84 

2.4 

37 

2.92 

2.0 

24.9 

3.68 

2.0 

21  .9 

4.63 

2.1 

20.2 

5.83 

1  .8 

15.4 

7.3 

1  .9 

14.8 

9.2 

2.6 

17.7 

11.6 

2.2 

13.5 

*  Rate  Constants  in  units  of  10"’  cm’ 
molec"'  s'\  and  are  accurate  to  ± 
0.5  X  10"’  cm'  molec"'  s" ' . 

••  Cross  sections  are  in  units  ofiO"’* 
cm^. 


Discussion 

These  results  can  now  be  Incorporated  into  a 
simple  model: 


0*  +  H.O  .f  H.O*  +0  ( 1  ) 

H.O*  +  H.O  -»  H,0*  +  OH  (3) 

(H.oL  +  hv  -*  H,0^  +  OH  +  e  (4) 
(HjO)j  +  0*  ->  H,0*  +  OH  ♦  0  (5) 

HjO*  4  Diffusion  (6) 

H,0*  ->  Diffusion  (7) 


Reactions  (1),  (3),  (6),  and  (7)  have  been 
proposed  and  discussed  in  previous  studies 


[Caledonia  et  al.|i967;  Hastings  et  al..  1988]. 
Reactions  (4)  and  (5)  have  not  been  suggested  as  a 
source  of  H,0*  in  the  shuttle  environment  previ¬ 
ously.  althoiuh  Reaction  (4)  has  been  postulated  a 
source  of  in  cometary  atmospheres  [Murad  and 

Bochsler.  1988].  Reaction  (3)  is  assumed  to  occur 
at  thermal  energy  at  a  temperature  of  •  273  K;  k. 
was  taken  to  be  2x10’^  cm’  molec"'  s*'  [Anicich 
and  Huntress.  1986]. 

As  discussed  in  another  study,  the  vaporiza¬ 
tion  of  water  generates  a  significant  quantity  of 
dimers  [Murad  and  Bochsler.  1987].  Observations 
made  on  the  shuttle  indicate  large  amounts  of 
particles  coincident  with  water  dumps  [Green  et 
al..  1987].  leading  to  the  euppoeitlon  that  the 
particles  are  ice  crystals.  Moreover,  meaeurements 
of  the  surface  temperature  of  the  shuttle  when  it 
is  in  darkness  indicate  it  to  vary  between  278°K 
and  238°K  [Humphries  and  Coyle.  1983].  From  these 
facts  we  can  assume  that  the  gaseous  water  envi¬ 
ronment  of  the  shuttle  contains  some  dimers,  since 
these  are  generated  in  the  equilibrium  vaporiza¬ 
tion  of  water;  in  this  analysis  we  take  the  dimer 
mole  fraction  to  be  0.12.  the  fraction  at  250°K 
[Murad  and  Bochsler,  1987].  The  photoionization  of 
the  water  dimer,  reaction  (4),  will  be  unimportant 
compared  with  either  (3)  or  (5).  since  its  esti¬ 
mated  rate,  is  6x10"^  s"'  [Murad  and  Bochsler, 
1987].  The  dominant  loss  term  for  H.O*  is  diffu¬ 
sion,  as  indicated  earlier  [Caledonia  et  al., 

1 987] . 

Since  k,  has  not  been  measured,  we  assume  that 
k^/k,  is  equal  to  the  ratio  of  the  calculated  rate 
constants  using  the  parametric  formula  of  Su  and 
Chesnavich  (1980): 

k  *  k,(0.3371P+0.62)  (7) 

with 

kj^«2»e(a/m^)'  and  P^g/CaAT)*  (8) 

where  k,  is  the  Langevin  rate  coefficient,  a  and  p 
are  the  polarizability  and  the  dipole  moment  of 
the  neutral  and  m  the  reactant  reduced  mass. 
Taking  p(H,0)  .  r.Sfp,  m( (H.0),.2.6D  and  ot((H,0),) 
«  <»(H,0)  »  1.54x10"  ®  cm’,  we  calculate  k  ■ 
3.06x10*’.  k  .  2.9x10*’,  and  k,  •  3.5x10^’  cm’ 
molec*'  8*',  leading  to  k  /k,  w  1.14.  At  a 
collision  energy  of  4.5  ev,  the  center-of-mass 
kinetic  energy  is  2.1  and  1 .4  eV  for  reactions  (1) 
and  (5),  respectively.  From  our  data  we  obtain 
k  *2.35x10*’  at  E_  •  2.1  eV  and  2.5x10*’  at  E_  * 
1.4  eV.  Then  at  eL  *  1.4  eV.  k  *  2.85x10*’  cm’ 
molec*'  s*',  sugge^ing  that  with  124  dimer  144  of 
the  H,0*  may  be  due  to  reaction  (5). 

Following  the  above  analysis  we  can  now  deduce 
the  steady  state  equations  for  [H^O*]  and  [H,0*]: 


k,[0*][Hj0]  *  k,[Hj0][H,0*]  +  (1/rj)[Hj,0*]  (10) 

and 

k,[H,0]  [H  O*  ]  +J  [ ( H  O ).]  ♦k  [0*  ]  [  ( HO)] 

’  *  *  *  (l/r,J[H,3*]  *  ®  (11) 


As  discussed  above,  the  second  term  in  equation 
(11),  the  photoionization  of  (H.O),,  can  be 
ignored,  and  we  take  [(HjO)^]*  o.  Ik  [H.O]. 
Equations  (10)  and  (11)  can  be  rewritten: 

[0*]/[Hj0*]  *  k,/k,  +  (1/rj)/k,[HjO]  =  t  (10’) 


and 


k,[H,0]  +  0.12  k,[0*][H  0]/[H  0*] 

’  '  ’  '  =(l/r,)[H,0‘]/[Hj0*] 


(11') 


Since  kj  fa  k,.  we  derive: 

[HjO]  a  (1/T^)/k,(«  -  1) 


and 

(l/rg)(k,+0.12«k,)/k,(«-1)  «  ( 1 /r,)[H,0*  ]/ [H,0*  ] 
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If  we  now  make  two  reasonable  assumptions,  namely 
kj  «  k,  and  «  r^,  we  obtain 

[H,0*]/[Hj0*]  =  (1  +  0.12«)/(«-l) 

Comparing  this  expression  with  the  ion  current 
ratio  of  0.1-0. 4  reported  by  Narcisi  et  al. 

[1984],  we  find  the  limiting  condition  S  <5 
corresponding  to  [H,0*  ]/ [H,0* ]  <  0.4.  There  is  no 
upper  limit  for  S  since  [ojO*  ]/ [HjO*  ]  >  0.12. 

The  limit  6  <  5  gives  an  upper  limit  for 
[HO],  depending  on  k  and  r..  With  r,  =0.1  s 
and  k,  =  2.5x10"’  mol  cm"’  s’  we  obtain  [H^O]  < 
ixio’  mol  cm"’  (-  3x10"®  torr)  which  is  somewhat 
lower  than  the  values  calculated  by  Caledonia  et 
al.  [1987]:  [HjO] =0 . 6-2 . 4x1 o’  mol  cm"’.  It  is, 
however,  consistent  with  the  values  calculated  by 
Wulf  and  von  Zahn  [1986]  from  neutral  mass  spec¬ 
trometer  measurements  ([HjO]  =  2  x  10®  cm"’  under 
quiescent  conditions  and  10'®  cm"’  during  RCS 
firings).  The  lower  value  of  [H^O]  calculated  in 
this  work  is  due  both  to  the  new  value  of  k,  and 
due  to  the  inclusion  of  reaction  (5).  the  reaction 
between  0*  and  water  dimers.  Of  course,  since 
water  concentration  is  inversely  proportional  to 
the  ion  drift  times,  and  r  .  shorter  drift 
times  would  lead  to  higher  water  pressures. 

We  conclude  from  this  study  that  the  partial 
pressure  of  water  is  likely  lower  than  previously 
thought,  suggesting  that  the  Space  Shuttle  envi¬ 
ronment  is  more  benign  for  optical  measurements 
than  previously  thought. 
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